Large-scale Dynamo Action Driven by Velocity Shear and Rotating Convection 
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By incorporating a large-scale shear flow into turbulent rotating convection, we show that a 
sufficiently strong shear can promote dynamo action in flows that in the absence of shear do not act 
as dynamos. Our results are consistent with a dynamo driven by either the shear-current effect or by 
the interaction between a fluctuating a-effect and the velocity shear; they are though inconsistent 
with either a classical a 2 or au mean field dynamo. 
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Magnetic fields are observed in virtually all cosmical 
bodies, from planets to stars and accretion discs; in many 
cases their presence can be categorically attributed to dy- 
namo action. The most pressing problem in astrophysical 
dynamo theory is to explain the generation of large-scale 
magnetic fields; i.e. fields with significant energy on scales 
large compared with those of the driving flow. The Sun, 
with its global magnetic field manifested through surface 
activity such as sunspots, represents the most well-known 
example of a large-scale dynamo. 

Astrophysical dynamos are often studied within the 
framework of mean field electrodynamics, a tremendously 
elegant theory that describes the evolution of a mean 
(large-scale) magnetic field in terms of transport coeffi- 
cients determined from averaged small-scale properties of 
the flow and field. The generation of magnetic field can 
then be ascribed to the a-effect, which relates the mean 
electromotive force to the mean magnetic field. The a- 
effect is non-zero only in flows that lack reflectional sym- 
metry [l[; consequently helical flows can be regarded as 
prime candidates for large-scale dynamo action. Indeed, 
in certain limiting cases the relation between a and helic- 
ity can be made explicit [HQ; however, and importantly 
there is no theory relating these two quantities when the 
magnetic Reynolds number Rm 3> 1 and the Strouhal 
number St « 1, the case of astrophysical relevance. Nu- 
merical simulations reveal that the relationship between 
a and helicity is indeed far from straightforward (3j. 

The most natural system for investigating astrophys- 
ical dynamo action is that of rotating thermal convec- 
tion [J, [5(. Recent studies of convection in a domain 
of large horizontal extent — namely one that encom- 
passes many convective cells — have demonstrated that 
although healthy dynamo action ensues provided that 
Rm is sufficiently large, there is no evidence of any signif- 
icant large-scale magnetic field in spite of the helical na- 
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ture of the convection @, 0] . Indeed, attempts to measure 
the a-effect directly reveal a strongly fluctuating quantity 
with a very small mean. The similarity of the spectra of 
the magnetic fields driven by rotating and non-rotating 
convection — for which the flows are not helical — pro- 
vides further evidence that the dynamo is controlled by 
small-scale processes (such as stretching and cancellation; 
see Q) and not by mean field processes (such as a lack 
of reflectional symmetry). 

The failure of rotating turbulent convection to act as 
a large-scale dynamo — as an a 2 dynamo in mean field 
parlance — suggests that the notion that helical flows 
will necessarily lead to large-scale field generation is too 
simplistic. It is, of course, the case though that most as- 
trophysical bodies possess a strong large-scale shear flow 
(differential rotation) and, indeed, most mean field as- 
trophysical dynamo models incorporate this feature. It 
is therefore of interest to examine the additional effects 
arising from incorporating such a shear into the rotat- 
ing convection model. It should though be stressed that 
these effects — if favourable in terms of dynamo action 
- should not disguise the inherent failure of the basic 
model to generate a large-scale field. 

One can envisage four possible beneficial effects of the 
shear on the mean field dynamo process: (i) that the 
large spatial scale of the shear leads to an enhanced a 
through greater spatial correlation of the small-scale mo- 
tions [7|, [9(; (ii) that even though the mean a remains 
small there may nonetheless be an effective aw dynamo 
when the shear is significant; (Hi) that the anisotropy 
induced by the shear may lead to a significant shear- 
current effect [13, [HI [H| i (} v ) that the shear may interact 
with temporal fluctuations in a to produce an effective 
mean field dynamo [H El]. In this paper we describe 
the effects of introducing large-scale velocity shear into 
the model of @, 0] in order to explore these various pos- 
sibilities. 

As in @, 0] we consider a plane Boussinesq convective 
layer (0 < x, y < A, < z < 1) with rotation about the 
vertical axis. In order to investigate turbulent dynamo 
action it is necessary to consider domains of large hor- 
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izontal extent, which encompass many convective cells; 
this though has to be balanced by the computational de- 
mands resulting from taking A large. Here we take A = 5, 
which results in 0(100) convective cells in the domain. 
This basic model is extended by the inclusion of a hori- 
zontal flow of the form 

U Q = C/ O cos (!) 

accomplished by replacing u with u+Uq in the governing 
equations. It should be noted that although a flow with a 
large-scale component (i.e. with the same spatial depen- 
dence as ([J)) does indeed result from this prescription, 
this component is not necessarily the 'target flow' given 
by CQ); the hydrodynamic state that ensues depends on 
interactions between the shear flow and convection and, 
possibly, on instabilities of the shear flow itself. It is also 
important to note that the scale of variation of this shear 
flow is much greater than all scales of the convection; this 
is essential if the results are to be explained within the 
mean field framework. Tobias et al. _15j have presented 
results in a related geometry but with a very different 
shear flow, namely one that has no horizontal structure 
but has a strong vertical variation. 

In order to elucidate the role of shear in this initial 
study we focus on the regime in which convection is fairly 
vigorous but in which there is no dynamo action in the 
absence of shear; specifically we set the Rayleigh number 
Ra = 150 000, the Taylor number Ta = 500 000, the 
Prandtl number = 1 and the magnetic Prandtl number = 
5; this leads to a Reynolds number 60 and a magnetic 
Reynolds number Rm « 300. A useful a priori measure 
of the imposed shear is given by the shear parameter S, 
defined by 

S = U (l/u rms L) , (2) 

where u rms is the rms velocity in the absence of shear, 
L is the scale of the shear and I is the horizontal scale 
of the convection cells in the absence of shear. For the 
parameters used here, S ~ C/o/300. It should be pointed 
out that there is also an effective value of S 1 , S e g say, 
defined analogously to ^ but involving the shear flow 
that emerges dynamically in the sheared covective state; 
this though can only be defined a posteriori. 

We have investigated flow and dynamo properties for 
the range < S < 7. A weak seed magnetic field of zero 
mean is introduced into an established, stationary, purely 
hydrodynamic state of sheared convection. Fig. [1] shows 
the evolution of the magnetic energy versus time for a 
range of values of S, and Fig. [2] the kinematic growth 
rate 7 as a function of S. We see immediately that dy- 
namo action ensues for sufficiently large values of S, al- 
though the dependence of 7 on S is not straightforward. 
Following the onset of dynamo action (with the critical 
value of S lying in the range l/3<5<l/2)7is linearly 
related to S, the strongest dependence possible (l6| . For 
larger S though this simple relationship no longer holds. 
This can be explained, at least partially, by inspection 
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FIG. 1: Magnetic energy evolution for a range of S. In terms 
of increasing linear growth rate, 5=1/3 (not a dynamo), 
2/3, 5/3, 5, 20/3, 10/3. 
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FIG. 2: Growth rates of the magnetic field versus S. 



of the purely hydrodynamic states. For the two largest 
values of S considered (S — 5, S — 20/3), there is a 
transition from the mode of convection that occurs for 
the smaller values of S; specifically, the proportion of en- 
ergy in the 'target mode' is much smaller, leading to a 
reduction in S c s- It is also interesting to note that for 
S > 1 the amplitude of the saturated magnetic energy is 
fairly insensitive to the value of S. 

Given that our motivation was the investigation of pos- 
sible large-scale dynamo action it is important to examine 
the spatial structure of the dynamo-generated magnetic 
fields, in both the kinematic and dynamic regimes. Fig. [3] 
shows the spectra of the horizontal fields for S — 5/3 and, 
for comparison, the spectra of the (small-scale) dynamo 
field for Ra = 1 000 000 in the absence of shear. Of par- 
ticular significance is that for the case of S = 5/3 there 
is roughly equal energy in all modes of scale comparable 
with and greater than that of the driving convective flow, 
in contrast to the case of no shear, for which the spec- 
trum is peaked at the scale of the convection. In both 
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FIG. 3: Horizontal power spectra for the magnetic field in 
both the kinematic (dashed) and dynamic (solid) regimes. In 
(a) S = 5/3, Ra = 150 000; in (6) S = 0, Ra = 1000 000; in 
both cases Ta — 500 000. The spectra were computed over the 
interior region of the domain (0.06 < z < 0.94). The arbitrary 
amplitudes of the kinematic spectra have been scaled so as to 
be on the same plot. 



cases we note the similarity of the shape of the spectra in 
the kinematic and dynamic regimes, indicating that the 
structure of the field in the final nonlinear state is deter- 
mined, to a large extent, from kinematic considerations. 

We have also directly determined the a-effect, by im- 
posing a uniform horizontal magnetic field and measur- 
ing the induced electromotive force. Since this procedure 
has an unambiguous interpretation only in the absence of 
small-scale dynamo action |17| . we have considered the 
value Uq = 100 (S ~ 1/3), which is strong enough to in- 
fluence the flow but is not quite strong enough to induce 
dynamo action. Fig. [Hz shows the time history of the lon- 
gitudinal a-effect (i.e. an calculated from £ x = auBo x ), 
obtained from a spatial average over half the domain [a] , 
for S = 0. As discussed in detail in Q, even though 
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FIG. 4: (a) Longitudinal a-effect versus time for S = 0; (b) 
a, the cumulative temporal average of a, for 5 = 0; (c) a for 
S = 1/3; (d) 57 for S = 1/3. 



a is the result of a spatial average over many convec- 
tive cells, it remains remarkably fluctuating in time, with 
large fluctuations and a small mean. Thus, as shown by 
the cumulative average in Fig. 3f>, a further long tem- 
poral average is needed in order to pin down a, with 
the resulting value being small in comparison with the 
rms velocity; from Fig. 2f) it can be seen that the long- 
time average value of a is given by a « 0.05, whereas 



4 



u rms w 60. For S = 1/3, when the influence of the shear 
on the flow is by no means negligible, there is essentially 
no difference in the behaviour of a to that when S = 0; 
as can be seen in Figs. Ht,d, it is again characterised by 
large fluctuations and the same small mean. 

Having thus shown that the introduction of a large- 
scale shear flow does indeed promote vigorous large-scale 
dynamo action, we should return to the four possibilities 
discussed earlier. From the considerations immediately 
above, revealing a to be essentially unchanged by the 
shear flow, we can rule out possibility (i). For a con- 
ventional auj dynamo model we might expect the growth 
rate to vary either as S 1 ' 2 (for a disturbance of fixed wave 
number) or S 2 / 3 (if the optimal wavenumber is permitted 
in the system) . Our calculations show that once dynamo 
action sets in then, for a range of S, the growth rate varies 
linearly with S; thus possibility (ii) is also not consistent 
with the results. Both remaining possibilities would seem 
to allow the growth rate to be linearly proportional to S 
[HI, Ulj]- However, distinguishing between the two is far 
from straightforward. Although the physical mechanisms 
are quite distinct — the shear-current effect depends on 
second order velocity correlations, whereas the fluctuat- 
ing a mechanism depends on fourth order correlations — 
they are both manifested as non-diffusive contributions 
to the turbulent diffusivity tensor fiijk- 

We have demonstrated conclusively that the flow re- 
sulting from the interaction of a large-scale shear flow 



and turbulent rotating convection can lead to large-scale 
dynamo action, i.e. the generation of magnetic fields with 
a significant component of energy on scales large com- 
pared with that of the convective cells. This may be 
significant in understanding the generation of large-scale 
fields in astrophysical bodies. We have concentrated here 
on the regime in which the convection, although fairly 
vigorous, does not induce dynamo action of itself. Thus 
the magnetic Reynolds numbers involved are fairly mod- 
est. Obviously it is also of interest to investigate the role 
of shear on the small-scale dynamo action that sets in 
at higher Rm and for which the underlying mechanism 
is not related to global measures, such as helicity, but 
instead to local stretching and folding properties of the 
flow, characterised, for example, by Lyapunov exponents 
and cancellation exponents [19j | . Our results on this will 
be presented in a future paper. 
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